RKKY interaction in SDW phase of iron-based superconductors 
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Using the multiband model we analyze the Ruderman-Kittel-Kasuya-Yosida (RKKY) interaction 
between the magnetic impurities in layered ferropnictide superconductors. In the normal state the 
interaction is spin isotropic and is dominated by the nesting features of the electron and hole bands 
separated by the antiferromagnetic momentum, Qy^^?. In the AF state the RKKY interaction maps 
into an effective anisotropic XXZ-type Heisenberg exchange model. The anisotropy originates from 
the breaking of the spin-rotational symmetry induced by the AF order and its strength depends 
on the size of the AF gap and the structure of the folded Fermi surface. We discuss our results in 
connection to the recent experiments. 

PACS numbers: 74.70.Xa, 75.30.Fv,75.30.Hx 



I. INTRODUCTION 

The oscillatory Ruderman-Kittel-Kasuya-Yosida 
(RKKY) interactioiy^ of the localized magnetic mo- 
ments in a metal has always played an important role 
in revealing the nature of the magnetic interaction 
in metals with partially unfilled d- and /-electron 
shells. This indirect interaction is the result of the 
spin polarization of conduction electrons produced 
by the exchange interaction of the localized moments 
with conduction electrons where the distance between 
two localized moments controls the strength of their 
effective exchange. Originally formulated for the three- 
dimensional spherical Fermi surface, RKKY interaction 
has been also analyzed for the two-dimensional (2D)^J^ 
as well as one-dimensional (1D)F electron gas. A partic- 
ular interesting situation arises in a highly anisotropic 
Fermi surface with nesting. In particular, it was shown 
that the RKKY interaction in such a case consists 
of several terms originating from flat regions in the 
fermionic spectrum (van Hove regions) and those which 
describe the interference between contributions from 
their vicinities.^ Note that the latter terms have an 
overall prefactor cos(Q • r) with Q = (7r,7r). In the 
nearly nested situation this last term is present down to 
the interatomic distances and favors the commensurate 
antiferromagnetic ordering of the localized moments. 

It is important to bear in mind that in order to evalu- 
ate the RKKY interaction at the distances r the details 
of the fermionic dispersion on a scale of 1/r in the /c-space 
is required. As a result the fine details of the Fermi sur- 
face are only important at largest distances while even 
the approximate knowledge of the electronic spectrum 
over the whole Brillouin zone is often enough for a good 
description of the RKKY interaction at the interatomic 
distances in r space. In this regard it is quite instructive 
to analyze the RKKY interaction in recently discovered 
Fe-based superconductors^. Band structure calculations^ 
and experimental probes such as angle-resolved photoe- 
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mission (ARPES) and quantum oscillatio: 
iments show that to a good approximation the Fermi 



surface topology of iron-based superconductors consists 
of the small sized circular hole pockets centered around 
the F— point (0,0), and elliptic electron pockets centered 
around the (±7r,0), and (0, ±7r)-points of the unfolded 
Brillouin zone (BZ). The pockets are nearly of the same 
size which results in the nesting properties of the elec- 
tron and hole bands at wave vectors, Qi ( Qi = (7r,0), 



and Q2 = (0,7r)), i.e. 



Given the elec- 



tronic structure of ferropnictides, it is natural to as- 
sume that magnetic order emerges, at least partly, due 
to near- nesting between the dispersions of holes and elec- 
tronPKni. 
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FIG. 1. (color online) Schematic Fermi surface topology of 
iron-based superconductors in the unfolded Brillouine zone 
(BZ) with 1 Fe per unit cell containing circular hole pocket 
centered around the F— point and two electron pockets cen- 
tered around the (7r,0) and (0,7r) points, respectively. 

Here we analyze the novel aspects of the RKKY in- 
teraction in iron-based superconductors which arise due 
to the peculiar Fermi surface topology in these systems. 
The origin of the local moments in ferropnictides can be 
either 4/-electrons in ReFeAsO series (Re- is a rare-earth 
element )P^^ or possibly partially localized d-electrons 
which arise due to proximity to a Mott insulatoi^^ which 
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FIG. 2. (color online) The contour mesh of the RKKY interaction for the normal (a) and the SDW state for Jx (b), and Jz 
interaction (c). In the SDW state we employ W — AQmeV ^ e = 0.5, and 5% of the electron doping. 



coexist with itinerant ones. Our primary interest is to 
investigate the evolution of the oscillatory behavior of 
the RKKY interaction in the presence of nesting in the 
normal and in the antiferromagnetic states of iron-based 
superconductors. Analyzing the RKKY interaction in 
the antiferromagnetic state we find spin space anisotropy 
of the interaction, a feature that has not been reported 
so far. We will also study in detail the influence of the 
model parameters like ellipticity of the electron pockets, 
and SDW gap size on the spatial variation of the RKKY 
interaction. 

The paper is organized as follows In Sec. II we eval- 
uate the RKKY interaction for a three band model and 
present its analytical form for SDW and normal state 
regimes. Using these results we evaluate the RKKY in- 
teraction numerically in Sec. Ill and discuss its relevance 
for the experiments. We finally present a summary and 
conclusion in Sec. IV. 



II. MULTI BAND RKKY INTERACTION IN 
THE NORMAL AND SPIN DENSITY WAVE 
STATES 

In this investigation we employ a minimal model of 
interacting 3d electrons and local moments in ferropnic- 
tides with a circular 3d hole Fermi surface (FS) centered 
around F-point (a-band) and two elliptical electron FS 
pockets centered around (±7r, 0) and (0, ±7r) points in 
the unfolded BZ (6-bands) (See Fig{l]). 

The Hamiltonian of the system of localized magnetic 
moment impurities in the multi band conduction electron 
sea is defined by 



H 



int •) 



(1) 



where 1-Lc = + conduction electron Hamil- 

tonian which is given by: 



k,Q!,cr 

k,<T 



.(2) 



Here, C^^,^ refers to the creation operators of the con- 
ducting electrons. In particular, a\^^ (^k7cr) creates an 
electron with spin cr in the hole (electron) band. The 
tight-binding energy dispersion of the electron and hole 
bands can be parametrized as follows 



£k = So + (cos kx + cos ky) 



(3) 



= ^0 + ([1 + e] cos{kx + tt) + [1 - e] cos{ky)) - fie 
^k = ^0 + ([1 - e] cos(/Ca;) + [1 + e] cos(/c^ + tt)) - /ie , 

where e accounts for the ellipticity of the electron pockets 
and £q = 0.05eV is the chemical potential for zero doping. 
Following our previous analysis^ we use the following 
hopping matrix elements th = 0.85 eV, te = —0.68 eV, 
/j^h = 1.44 eV, and /ig = —1.23 eV which accounts for 
the Fermi velocities and sizes of the Fermi pockets, see 
Ref.E) 

The interaction part of the conduction electron Hamil- 
tonian contains density-density interactions between hole 
and electron bands which give rise to a SDW order be- 
tween the hole pocket and one of the electron pocket 
located around the (tt, 0) point of the BZ^^. Assuming 
the experimentally observed Qaf — Qi = ('^r, 0) SDW 
ordering wave vector a standard mean-field decoupling 
yields the self-consistency condition for the SDW order 
parameter W oc Xlk(^k5^ik5'^55'). The resulting mean- 
field Hamiltonian has the form 



k<7 



H.c. 



, (4) 



where the spin index a = ±1 refers to the spin \ and | 
respectively. Now applying the unitary transformations 



ttka = ^kCka - ^k<^ka 
^Ik+Qia = Cr[^kCka + ^k^^ka] 



(5) 



the mean-field Hamiltonain for the conducting electrons 
can be diagonalized and the coefficients of the transfor- 
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mation are given by 
. 1 



(6) 



The diagonalized Hamiltonian has the form 



k,cr 



(7) 
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Jn 
Jx 



where J ex is the exchange couphng constant. It can be 
obtained from a more microscopic Anderson type model 
introduced in Ref. |24j Here is the moment of localized 
/-electrons at site R^, and s(r) is the spin of conduction 
electrons. 

In the spin density wave state we employ the standard 
second order perturbation theory with respect to l-Lint- 
Its application is straightforward and one finds after some 
algebra the RKKY interaction which describes the inter- 
action between two local impurity spins at the positions 
i and j in the form of an XXZ type effective exchange 
Hamiltonian 
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FIG. 5. (color online) Oscillatory behavior of the RKKY 
interaction along x direction in the normal (Jn) and SDW 
[Jx and Jz) states. 



^RKKY — Jx i^i^j + ^i^j) + Jl^ S^Sj (10) 

The magnetic anisotropy of the form Jx = Jy 7^ Jz iii 
this expression appears through the SDW order which 
is polarized along z-spin quantization axis. Specifying 
Ri = and Rj = r = {x^y) these effective exchange 
couplings are given by 
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FIG. 6. (color online) Jx and Jz components of the XXZ 
RKKY interaction in the SDW state along two crystallo- 
graphic directions. As above we employ W = 40meV, e = 0.5, 
and 5% e-doping. 
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Note that the electron pocket located around the (0, tt) 
point of the BZ remains intact and is not involved in 
the SDW formation. In addition, the SDW order explic- 
itly breaks the spin-rotational symmetry. As a result the 
RKKY interaction will not be spin-rotationally symmet- 
ric and will have a distinct magnetic anisotropy induced 
in the SDW state. 

The interaction part of the local magnetic moments 
with conduction electrons is given by a contact exchange 
term 
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Here /(k, 7) is the Fermi function, and the SDW coher- 
ence factors, ^kWcr'' defined as 
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Setting = it is easy to verify that in the paramag- 
netic or normal state regime the RKKY interaction sim- 
plifies to the usual expression in two-dimensional metals 
with nesting properties^ 



RKKY 
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FIG. 7. (color online) The effect of the ellipticity on the modified XXZ Heisenberg type RKKY interaction in the SDW state, 
Jcc, along (100) (a) - and (010) (b) - crystallographic direction, respectively. 



where the interaction is now isotropic in the spin space 
{Jx = Jy = Jz = Jn)- The effective exchange couplings 
are then given by 

77' \77' 

(15) 

here Qe^ = Qi, Qes = Q2, Q/i = 0, and x^^'(r) is mag- 
netic spin susceptibility of conduction electrons (Lind- 
hard response function) which is given by 



^77' (r) = ^e^(k-k')-r 



/(k,7)-/(k',Y) 

^k ^k' 



. (16) 



III. NUMERICAL RESULTS AND DISCUSSION 

Based on the equations above we present in the follow- 
ing the results of the numerical evaluation of the RKKY 
interaction in the SDW and the normal state phases. We 
first present an overall behavior of the RKKY interac- 
tion. In particular, Fig. ([2]) shows the contour mesh of the 
RKKY interaction for 5% electron doping and ellipticity 
parameter e = 0.5 as a function of interatomic distances 
for both SDW and normal state regimes. In addition to 
the breaking of the spin rotational symmetry, there is an- 
other important difference in the behavior of the RKKY 
interaction between the normal and the SDW state of 
iron-based superconductors. This concerns the absence 
of tetragonal symmetry in the SDW state. As clearly 
seen by comparing Fig.([2|(a) and Fig.([2|(b)-(c) the C4 
symmetry present in the normal state is broken down to 
C2 symmetry. Apart from these differences the RKKY 
interaction show also some similarities. In particular, 
close to the impurity position the effective interaction is 



ferromagnetic (FM) in both cases and becomes antifer- 
romagnetic(AF) at a distance comparable to the lattice 
constant. In the asymptotic regime {x/a^y/a ^ 1) the 
oscillatory behavior from FM to AF and vice versa sets 
in. As expected the amplitude of oscillations decreases 
with increasing distance between the local moment im- 
purities. 

For clarification of the role played by the interband and 
intraband scattering we display in Fig. [S] the contribution 
of each term separately in the normal state. These fig- 
ures show that the inter band contribution to the RKKY 
interaction is much larger at shorter distances than the 
intraband one. This arises again due to pronounced nest- 
ing features of the electron and hole bands. In addition, 
observe that a finite ellipticity introduces some asymme- 
try along X and y direction for each component of the 
RKKY interaction where electron pockets are involved. 
This asymmetry, however, averages out in the full RKKY 
interaction which possesses again the tetragonal symme- 
try. 

This is not any longer the case in the SDW state 
as shown in Figj4] for Jx component. Here, the intra- 
and inter-band contributions to the RKKY are strongly 
anisotropic. This originates from the fact that the SDW 
has Qi = (tt, 0) ordering wave vector and as a result of 
its ordering the tetragonal symmetry is broken. It is in- 
teresting to notice that such effect was found recently in 
EuFe2As2 where the magnetic anisotropy of the Eu mag- 
netic moments was changing across the SDW transition 
temperature^^. This change in the magnetic anisotropy 
is in direct agreement with our results. Note also that 
the sum of all contributions is anisotropic along a and 
h crystallographic directions in both Jx and Jz compo- 
nents of the interactions. In particular, as clearly shown 
in Fig 2|b)-(c) the ferromagnetic interaction is extended 
in both Jx and Jz along the y direction which is perpen- 
dicular to the AF ordering of the Fe-plane. 

To see the changes on the quantitative level we show 
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FIG. 8. (color online) Influence of the SDW gap magnitude on the RKKY interaction parameter, J^, for (lOO)-direction (a); 
and for (OlO)-direction (b). 



in Fig|5]the dependence of the RKKY interaction along 
X direction in both normal and SDW states. We ob- 
served that the amplitude of the oscillations increases in 
the SDW state but overall the dependence remains the 
same. Namely, it is ferromagnetic for short distances 
and then oscillates between positive (AF) and negative 
(FM) values. In addition, Jz and Jx RKKY interac- 
tions have slightly different period which results in the 
fact that they may have opposite signs for a given inter- 
impurity distance. This difference is associated with the 
structure of the SDW matrix elements which appear in 
the and components of the RKKY interaction dif- 
ferently (see Eqs.(ll)-(12)). An additional effect of the 
SDW is shown in Figj6] where we plot the behavior of 
Jz and Jx along a and h crystallographic directions. As 
clearly seen from this figure, the period of the oscillation 
is not only different for Jx and Jz components but also 
for each of them in the x and y crystallographic direc- 
tions. In particular, the antiferromagnetic XXZ regime 
with negative anisotropy at x :^ 2a is observed along 
(0,1) direction. On the other hand for (l,0)-direction 
the XXZ ferromagnetic Heisenberg model is dominant 
and the effective interaction changes to AF behavior only 
around y 2.5a. Thus through SDW spin-space and 
real space anisotropics are correlated. Furthermore along 
(l,0)-direction the magnetic anisotropy changes sign for 
the first time at a distance of about y 0.75a and it pre- 
vails for a longer period as compared to (0,1) direction. 
The nature of this difference is both SDW order and the 
structure of the remaining small pockets that occur due 
to folding of hole and one electron pocket located at (tt, 0) 
point of the BZ. Due to larger hopping of this electron 
band along x-direction the values of the folded bands 
in the SDW state are unequal along x and ^/-direction 
which is then reflected in the periodicity of the RKKY 
interaction. 

Furthermore we show the effect of the ellipticity in 
Fig j7| where one could clearly observe the increasing pe- 



riod of the oscillations for larger values of e. The same 
effect is observed for increasing electron doping. This is 
natural as both ellipticity and doping for a given value 
of W make the remnant pockets and the corresponding 
values of kp along x and y direction larger. This effect 
is almost absent in the normal state (not shown) which 
points out that in the SDW state the interaction between 
magnetic impurities will be strongly modified. 

The natural question arises whether the effects of the 
SDW state on the RKKY interaction become more pro- 
nounced for increasing size of the SDW gap and cor- 
responding increase of the magnetic moment. In FigjS] 
we show the evolution of the oscillatory behavior of the 
RKKY Jx interaction for different value of W. Note that 
the amplitude of the oscillations weakens upon increase of 
the SDW gap. This is due to the shrinkage of the remnant 
electron and hole FS pockets which arise due to folding 
of the BZ in the SDW. The larger becomes the SDW gap 
the smaller will be the remnant pocket size. As a matter 
of fact for some critical value of W the pockets involved 
in the SDW completely disappear from the Fermi surface. 
Therefore the only contribution to the RKKY will arise 
in this case due to electron pocket located at (0,7r), not 
involved in the SDW formation. This explains why Jz 
shows weaker oscillations along x-direction, while along 
^/-direction the oscillations are almost the same. This is 
because the oscillatory behavior of Jz along (01) direc- 
tion is determined by the electron pocket which remains 
intact in the SDW state while along (1,0) direction the 
SDW order gaps completely the FS and only slight oscil- 
lations are still visible. 



IV. SUMMARY AND CONCLUSION 

In conclusion we analyze the changes of the RKKY 
interaction in the SDW state of iron-based superconduc- 
tors. The generalized RKKY interaction in these com- 
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pounds is of an effective XXZ Heisenberg-type where the 
0(3) symmetry is broken but U{1) symmetry of the in- 
teraction for rotation around an axis parahel to the SDW 
polarization vector Wz is stih preserved. We show that 
for smah distances between the local moments, r < a, 
the interaction between local spins is ferromagnetic but 
for larger r it oscillates between AF and FM regimes with 
different periods and amplitudes which depend strongly 
on ellipticity or doping of electron pockets. In addition, 
the period of the oscillation strongly depends on the mag- 
nitude of the SDW gap and the structure of the Fermi 
surface in the folded BZ. 

Our main observation is that the RKKY interaction 
between magnetic impurities in SDW state become 
anisotropic below Tsdw- As a result, the magnetization 
of the rare-earth magnetic moments, already anisotropic 
by itself due to crystalline electric field effects, will 
experience additional temperature dependent anisotropy 
induced by the conduction electrons below Tsdw- Quite 
generally the effect of SDW ordering of Fe spins o n the 
rare-earth subsystem was found in several studieP^^^. 
However, the effect of induced anisotropy below Tsdw 
on the rare-earth magnetization was observed only 



recently in EuFe2(Asi_a^Pa^)2 system by measuring 
magnetic anisotropy of the Eu^+ ions above and below 
Tsdw- Iu particular, it was found that upon decreasing 
temperature the ratio the magnetization anisotropy of 
Eu spins, Mab/Mc, becomes temperature dependent 
below Tsdw refiecting the infiuence of the SDW oide^^. 
This is in direct agreement with our results. Our further 
observation that the magnetic anisotropy is then also 
reflected in the spatial anisotropy was not yet observed 
as it requires the use of untwinned crystals or the use 
of the local probes such as nuclear magnetic resonance 
(NMR). It would be interesting to check this effect 
experimentally. 
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